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INTRODUCTION 

The Grant NGR-06-002-147, "Theoretical Investigations of Plasma 
Processes In the Ion Boobardment Thruster," Is concerned with 1) the 
sputtering of the accelerating grid, and the cathodes of the hollow 
cathode and neutraliser discharges, ii) the deposition of the sputtered 
atoms on system components such as the solar energy collectors, and 
iii) the analysis of the thruster discharge. The progress made on 
these subject in the period from 7.1.74 to 8.30.75 is communicated herein 
A physical model for a thruster discharge is developed consisting 
of spatially diverging plasma sustained electrically between a small 
ring cathode and a larger ring anode in a cylindrical chamber with an 
axial magnetic field. The associated boundary-value problem for the 
coupled partial differential equations with mixed boundary conditions 
is solved in closed form for the electric potential, the electric 
field, current density, and velocity distribution. As a result of the 
Lorentz forces, the plasma rotates with speeds of the order 10 cm/sec 
around the chamber axis. It is shown that at sufficiently large Hall 
coefficients and/or Hartmann numbers i) the radial spreading of the 
discharge is reduced by the external magnetic field, and ii) the dis- 
charge fields exhibit an oscillatory spatial structure off the discharge 
axis, e.g. , current flows in alternating axial directions. 

By means of quantum-mechanical perturbation theory, a formula for 
the number S(E) of atoms sputtered on the average by an ion of 
energy E is derived from first principles. The theory agrees with 
experimental sputtering data in the low energy region above the thres- 
hold, and leads to the correct sputtering thresholds. As an appli- 
cation mercury-metal atom scattering cross sections are determined 


I 


! 
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by quantitative comparison of the theortetlcal and experimental S(E)-values 
for sputtering mercury Ions and polycrystalline target materials, such 
as Ag, Au, Co. Cu, Fe, Mo.Nb.Pt, Ta. Ti, W, and Zr. 

The boundary^value problem describing the diffusion of the sputtered 
atoms through the surrounding rarefied electron-ion plasma of ion pro- 
pulsion systems to the system surfaces is formulated and treated 
analytically. It is shown that outer boundary-value problems of this 
type lead to a complicated Integral equation. Numerical results can 
be obtained by a considerable computer effort. 

The investigations reported herein represent preliminary 
communications. The final version of this work will be communicated 
in form of publications at a later date. 
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THEORY OF ROTATING DISCHARGE IN MAGNETIC FIELD 


ABSTRACT 

A physical model for a thruster discharge Is developed consisting 
of a spatially diverging plasma sustained electrically between a small 
ring cathode and a larger ring anode in a cylindrical chamber with an 
axial magnetic field. The associated boundary-value problem for the 
coupled partial differential equations with mixed boundary conditions, 
which describe the electric potential and the plasma velocity fields, 
is solved in closed form. The electric field, current density, and 
velocity distributions are discussed in terms of the Hartmann number H 
and the Hall coefficient ut. It is shown that the plasma fields 
exhibit an oscillatory radial structure at sufficiently large magnetic 
Interaction parameters H and «t. As a result of the Lorent* forces, 
the plasma rotates with speeds as high as 10^ cm/sec around its axis of 
symmetry at typical conditions. 


CATHODE 

I 

-2 

FIG. 0: Geometry of Spatially Diverging Discharge Between 

Cathode (Rj^) and Anode (R 2 ) in a Cylinder (R^) 
with Axial Magnetic Field (R-^R, ) . 
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In the ion thruster, the plasoa is produced by means of a diffuse, 
spatially diverging electrical discharge between a hollow cathode of 
small radius and a circular anode of considerably larger radius R2 

(at the cylindrical chamber wall).^^ The field lines of the electric 
current density J and of the confining magnetic field S cross under 
a nonvanishing angle (except at the chamber axis) so that the resultant 
Lorentz force 3 ^ x $ rotates the discharge around its axis of synmetry. 

In steady state, the magnetic body forces in azimuthal directions are 
balanced by viscous forces (boundary layers at the chamber walls) . 
Schematically, this is illustrated in Fig. 0 for a much simpler model 
of a discharge between a cathode (Ri) and an anode (R2) in a homogeneous 
magnetic field 5 ^ (R2 » Ri , rotation in direction u) . 

In view of the complicated geometry of actual thruster discharges 

1 2 ) 

and their inhomogeneous magnetic confinement fields, * the steady 
state rotation of a spatially diverging discharge in an external 
magnetic field is analyzed by means of the electrical discharge model 
in a homogeneous axial magnetic field depicted in Fig. 0 . The 

analysis is based on the magnetogasdynamic approximation, in which two 
characteristic nondisienslonal parameters occur, the Hartman number H 
and the Hall coefficient ut, 

00 o 

The symbols Introduced designate the electrical conductivity (0), the 
viscosity (m), the electron gyration frequency (w), and the electron 
moMntum relaxation tisie (t). H and oit are a measure for the strength 
of the Lorentz force relative to the viscous friction force and for the 
reduction of the current flow transverse to the magnetic field ^ , 


/ 1 I I 1 i I I I 

6 

respectively. Very little Is known about the transport coefficients 

0, y, and t for low pressure discharges, except the qualitative 

experiaental results that o Is anotaalously mall, y Is anomalously 

3) 

large, and t Is anomalously small. 

For the above reasons, the dynamics of the rotating discharge In 
an axial magnetic field will be discussed In terms of the nondlmenslonal 
parameters H and ut which are treated as variables within their 
value domain of practical Interest. H 2 1 and wt ^ 1. H and ut 
represent phenomenological parameters since a, y, and t have to be 
obtained by measuremnts. 
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THEORETICAL FORMULATIOM 


For a purely azimuthal flow field, v ■ {0, v(r,z),0}, the plasma 

behaves Incompressible, V • v ■ 0. From the continuity equation, 

V • (pv) ■ V • 7p ■ 0, it follows then that the density gradient 7p 

is everywhere perpendicular to the flow field v. These simple conditions 

are only approximately realized since the intensity of the secondary 

— — 4 ) 

flows Increases with increasing Reynolds number R • p„vR^/p. In 

n o 

absence of secondary flows, momentum cannot be balanced completely in 
the z-dlrectlon, since 

3p/&z ■ 0, for V - {0, v(r,z),0}, 

in accordance with the z -component of the magnetogasdynamic equation of 

motion. A strictly z-independent pressure field is physically not 

possible, since the axial pressure gradient across the boundary layers 

at the end plates z ±c is nonzero, although generally small, 

3p/3z A 0 (boundary-layer approximation for "transverse** pressure 
4) 

gradient) . 

In accordance with the magnetogasdynamic equations,^^ Ohm's law 
with Hall effect,^^ and the conservation equation for the electric 
charge density (V * 3^ * 0) , the rotating discharge in a hosu>geneous 
magnetic field is described by the following boundary-value 
problem for the azimuthal velocity v(r,z) and electric potential 


^(r,z) fields (Fig. 0; secondary flows neglected): 

4 " • If ®o^' If ^ ®o> 



( 1 ) 

( 2 ) 

O) 


f 


1 


I 


1 
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where 


v(r,z)^_j^ - 0 . -c < z < +c 
and 


(A) 

(5) 


-o [3^(r,z)/3z3^^j^ ■ I 6(r - Ri,2)^2irr , (6) 

[3^(r,z)/3r]^_j^ • 0 , -c < z < +c . (7) 

o 

The boxmdary conditions (4), (5), and (7) consider that the plasaia 

does not slip at the walls r ■ R and z ■ ±c , and that no current 
^ o 

flows into the cylinder wall r ■ R^, respectively. The bounlary 
conditions in Eq. (6) imply that the cathode (R;) and anode (R 2 ) are 
ring electrodes of vanishing radial width, Ar -► 0[6(r - Ri^^^/^wr ■ 
radial Dirac function] . The net current flowing through the discharge 
is by Eq. (6) 

R© \ 

-2tto f ■ * rdr - I J 6(r-Ri 2 ) ^ ^ 0 * 

0 0 * 


since the positive current (1 < 0) flows from the anode to the cathode 
(Fig. 0). In Eqs. (1) - (3), the transverse conductivity is given by 
Oj^ • 0/(1 •¥ . ( 8 ) 

Owing to the disregard of secondary flows, the boundary>value 
problem is linear since v(r,z) and ^(r,z) are described by the 
linear equations (2) - (7) . The nonlinear aquation (1) determines the 
pressure field p(r,z) which does not occur in Eqs. (2) “ (3). It 
should be noted that the induced magnetic field has been disregarded 
under the assumption that the magnetic Reynolds nusd>er is small. 


R « 1 


1 


,.4. 
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ANALYTICAL SOLUTION 


The cheraccerlstlc nondiaenelonel peraaetcre of the aegnetogaedynealc 
dlacherge problea under consideration are obtained by introducing the 
nondlMnslonel Independent end dependent variables. 


p • r/R , 0 < p < 1 • (9) 

0 “ 

Q ■ */c » ”1 i C £ +1 * 

and 

V(p,C) ■ v(r,s)/v^, ♦(p,C) - ♦(r.z)/*^ ♦ (11) 

where 

V 5 ♦ /R B , ^ 5 Ic/2voR2 . (12) 

o o o o o O 

In teras of the nondlaenslonal space variables and fields, the 
boundary*value problea deflried In Eqs* (2) - (7) assusMS for V(p,() 
and *(p,C) the fora: 





(13) 




-h2 ^ 


(14) 


where 

'^<P*^>p-l ■ 0 , -1 i C < +1 » (15) 

V(p,C)^,*l - 0 , 0<p<l . (16) 

and 

-[a6(p.C)/3Cl^.jl ■ «(P-P2.l)/P * 

l3A(p,C)/3plp.i - 0 , (18) 

with P 2 ^i = R 2 ,i/R^. The nondlaenslonal constants M, N, and are 

defined by 


I 


I 
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M-2 - (1 + w2t2)(R^/c)2 , N"2 - <Rj,/c)2 , 

**1 " U^T^) 


(19) 

( 20 ) 


In view of the similarity of the left sides of Eqs, (13) - (14) 
with Bessel's differential equation, + P“^Z^ + (k^ - p~^m^)Zj^ ■ 0, 
for cylinder functions Z (k p), partial solutions of the coupled 
Inhomogeneous equations are sought in the form. 


♦^(P.O - J^(k^p)f(0 , 


V^(p,c) - Ji(k^p)g(;) , 


( 21 ) 

(22) 


where Jo(k^p) - -Ji(k^p) and J[(k^p) + (k^p)“lji(k^p) • Jj,(V^* 
Substitution of Eqs. (21) - (22) into Eqs. (13) - (14) yields 

d2f 


dc" ■ “5 'v ■ \ »v 
^ . (kj + - k^ . 


(23) 


(24) 


where the eigen-values k^ > 0 are determined by the boundary conditions 
(15) ai^ (18) as the real roots of the transcendental equation. 


Jj(k^) "0, V* 1,2,3,.... 

Thus, the general solution of the coupled equations (13) 
by linear superposition as: 

m 

♦ (P.c) - I . 

V»1 


(25) 

(14) obtains 

(26) 


V(p,C) ■ I Ji(k p)g (O 


( 27 ) 


v-1 


By decoupling Eqs. (23) - (24) one finds for f^(C) 
differential equations of 4th order. 
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d^f 


d2f 


d?r - + n2) + n2h|1 5^+ - 0 


d*»g 


d^g 


with 


5^ - (kJ(M2 + n2) + + *^ . 0 






( 2 ») 

(29) 

(30) 

(31) 


as brundary conditions by Eqs. (16) - (17). In deriving Eq. (30), the 
Dirac function in Eq. (17) has been expanded. 


«(p - P2,l)/P ■ 2 J (\P) • 


(32) 


v-1 


In addition to Eqs (28) - (31), f^(C) And 8^(C) have to satisfy also 


the uncoupled Eqs. (23) - (24). With 


U) =U)..U,. i(t) .U) = '*4d . (tl s ~iH I 

IV V+* 2 V V-* 3 V V+* 4 V V-* 


(33) 


“v± * UkJ(M2 + n2) + n2h2]2 

- 41cJm2n2 )**}]** , (34) 

the general solutions for f^(C) ' e“^ and g^(C) •’ e“^ of Eqs. (28) • 
(29), can be written as: 

cosh 


slnh u. C 

f (O ■ A, — + B, . 

V Iv sinh ( 1 ), iv cosh u, 

Iv Iv 


slnh cosh 

^ ^2v slnh ^ ®2v cosh w_ * 
2v 2v 


sinh w. c cosh u. c 

g^(C) “ cosh 

sinh Wj c cosh w- c 

^ C ^ p 

2v sinh b> 2 ^ 2v cosh 


(35) 


( 36 ) 
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Only four of the eight Integration constants for any 

V ^ 1 are Independent; by Eqs (23) - (24), 

(“I - 

(“L - ■ ’'v‘*"‘=2v 

(37) 

(u? - k2M2)B, , - k M^D, 

' Iv V Iv V Iv 

(o)2 - k2M^)B, - k M^D, 

2V V 2v V 2v 

and 

[u»2 - (k2 + H?)n 2]C, - k n2h?A, 

■■ Iv V ^ Iv V Iv 

[«2 - (k2 + H?)n 2]C, - k n2h?A- 

2v V ' 2v V 2v 

9 ^ (3®) 

l“!v - 

I“lv - <“5 + ■ “/"“pav 


where the coefficient determinants of the pairs of corresponding 
equations in Eqs. (37) and (38) vanish owing to Eqs. (33) (34). 

Upon application of the four relations in Eq. (38), which are 


equivalent: to Eq. (37) by Eqs. (33) ^ (34), and the boundary conditions 
(31) , which give 


-C « C. = C , -D.. - D ED 

2v Iv V 2v iv V 


Equations (35) > (36) become: 


^v^^^ " k^N^H^ j”lv slnh " ^2v slnh w 


sinh 


slnh o» 2 ^C 


V 

iT5^ 


cosh w c 

0, - - 0 


Iv cosh 


2v f 
cosh w C 


u 


Iv 


2v cosh u, 


2v 


(39) 


(40) 
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where 


[ sinh «j c ainh w c I 
einh “ sinh 

[ cosh b), C cosh (D. C 1 

— r^\ . 

cosh u). cosh (i >2 J 


“iv = “iv - ‘“5 + 1 - 1-2 • 


( 41 ) 


(42) 


The boundary conditions (30) applied to Eq. (40) yield 

[J (k p.) + J (k p«)l 


k N^H? 

I 


j2(k ) [(I), n, cthu, -w^ 0_ cth(i>^ ) 

o V Iv Iv Iv 2v 2v 2v' 


(43) 


k N^H| 

L_ 

J2(k ) 

O V 


l“lv“iv'*'“iv-“2v“2v'*^2v> 


(44) 


Substitution of Eqs. (43) - (44) into Eqs. (40) - (41) gives as 
solutions for f^(C) and g^(C) • 


f (C) • J2(k ) 

V O V 


■ “2v“2v"*"2v’ 


[“ 

[• 


sinh (t), C 
Iv 


sinh (I),. c 


Iv sinh <i> 


- n 


Iv 


cosh 


2v sinh 


^1 

"2V J ‘“lv®l 


‘•'o<Vl> - •’o<V2” 


Iv cosh w 


- 0 


iv 


cosh u> 
2V cosh 


, 0, tglui). -( 1 )^ 0. tghu. ] 
Iv Iv • Iv 2v 2v ® 2v 

"2V J 


and 
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g„U) 




l“w«lv‘''"lv-“ 2 v“lv“'“lvJ 


[ 


'ainh Wj^C 
“iv 

ly^ 

Iv 


8inh W-. C 


slnh 
cosh 


'2v' 1 * tJptk'.p,) - J„(Vi)l 
“ 2 vJ l“lv“l/»'“lV-“2vV*'“2v' 


fcosh cosh “| 

[ cosh <i)j cosh a> 2 ^ J 


(46) 


Equations (45) - (46) form, together with Eqs. (26) - (27), the closed 
form solution of the problem of the rotating gai discharge in an axial 
magnetic field 

. J_^(k^p) ( Uo(k^Pl) + Jo<V2>i _ 


■ ■ vii \ 1“iv“iv«'“ip-“2Av“'“2v> 

sinh u„ 


o V 


[■ 


sinh 0 )^^; 


Iv sinh u 


- Si 


'2v sinh 


“2V^ ~ 
“2m . 


IV ■ 2v 


Foi Si tehb) ~u Si_ tghu)^,, 
^ iv iv * IV 2V 2v * 2v 




cosh 


cosh u„. t 


Iv cosh u> 


- Si 


iv 


2v 


cosh 


2V 

2v J 


2v 

(47) 


and 


- V(p,0 -N^Hf I " 




ivil J2(k) b“lV«lv"^^lV-“2vSv"^^“2vl 


[■ 


o' V 


sinh sinh 


sinh u) 


sinh 


^1 

J 


Iv 2V 

lJp<Vi> ' •’p<Wl 


rcosh 

cosh I 

[cosh 

" cosh J j 


(49) 


■ [“lv“lv‘»*“lV - “2v®2v‘»'“2vi 

Th» rpuinlng nondl»enplon»l ditehprgt fltlds 9 - - 

J - " 5/3 are given in terms of the solutions for ♦(p,C) snd V(p,;): 
o 

B* • -N 36/ap , E* ■ 0 , E* - -a*/c 


9 (49) 
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N(-|i+V), J, 


1 + 


U)T 

UTT^ 




li 

3 ; 

(50) 


where ** " 

If the cathode la in the plane a - -c (c • -D the anode is 

in the plane z - +c (C - +D. then the reference fields and 

[Eq. (12)] are negative, since I<0. The results are also applicable 

to the case where the anode is in the plane z - -c (C • -1) and the 

cathode is in the plane z - +c (C - +D. m the latter situation, the 

reference fields v and [Eq. (12)] are positive, since I>0. 

These explanations hold for magnetic fields pointing in the positive 

z- direction. B >0; v changes its sign if B <0 [Eq. (12)]. 

0 0 ^ 
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NUMERICAL ILLUSTRATIONS 

As an illustration, the radial (p) dependence of the nondlmenslonal 
discharge fields V(p,0. *(P.C). E^(P»0, E^(p,C). and J^(p,C) has 
been calculated for I < 0 in the cross sectional planes C - -0.9 
(cathode region), C ■ 0 (central region), and C ■ +0.9 (anode region) 
based on Eqs. (47)-(50). The remaining fields Jg(p,C) and J^(p.C) 
are simply proportional to J^(p»C) and E^(p,C)» respectively 
[Eq. (50)]. The characteristic (nondlmenslonal) magnetic Interaction 
numbers are treated as parameters: 

WT - 1, 10, 100; H - 1, 10, 100. 

-2 2 2 

The geometry parameter N is taken to be N • 1 so that M ■ 1 +w t , 
corresponding to R^ - c [Eq. (20)]. The radial positions of the 
cathode and anode are assumed to be; 

p^ - 0.01 (Rj^ - 0.01 R^) ; P2 - 0-9 (^2 “ V * 

The dimensional fields are negative everywhere where the nondlmenslonal 

fields are positive, and vice-versa [Eq. (11)] since v^ < 0 and 

a <0 for I < 0 [Eq. (12)]. 
o 

Central Region, c ■ 0; In Figs. 1-9, the potential ♦, the radial 
electric field E^, and the aslmuthal flow field V are represented 
versus p for C • 0 and ux - 1, 10, 100; H - 1, 10, 100. The 
extrema of V and E^* move towards the center p ■ 0 as wx and/or 
H are increased. ♦ always decreases from a maximum at P - 0 to a minimum at 
p ■ 1 [E (p • If C) • 0]. It is remarkable that at certain values of 
«x and H the fields ♦, E^ and V exhibit an "oscillatory" radUl 
distribution (E^ <0, V <. 0; the oscillation amplitudes of ♦ are 
generally too small to be visible in the ♦ curves) . 
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The Figs. 10-12 show the axial electric field versus p 

for t - 0 and wt - 1, 10, 100; H - 1, 10, 100. It is seen that 

is particularly strong in the center, p ■ 0, and changes repeatedly its 

sign at larger p values, e.g. in an oscillatory manner for larger 

uT and/or H values. This means that the bulk axial current flows in 

a narrow region, p * 0, in the direction anode-*«athode, outside of 

which the axial current density J, changes periodically its sign. The 

pronounced maximum of E at p * 0.9, is due to the effect of the 

z 

anode at p - 0.9, C - +1 on E^ in the plane C - 0. 

The Figs. 13-15 represent the radial current density versus 

p for t - 0 and ut - 1, 10, 100; K » 1, 10, 100. vanishes 

always for p • 0 (symmetry) and p ■ 1 (J^(p • 1» C) ■ 0] . The bulk 

of the radial current flows in a restricted radial region off the center 

p ■ 0 (c * 0!). The extremum of at p * 0.85 reflects the 

influence of the anode at p - 0.9, C “ +1 on in the plane C • 0. 

At sufficiently large values of ut and/or H, flows periodically 

forward and back radially in accordance with the oscillatory radial 

structures of E^ and the induced electric field - V. 

Cathode Region, e ■ -0.9; In Figs. 16-24, the fields ♦, E^, 

and V are shown versus p for C ■ -0.9 and ut ■ 1, 10, 100; 

H - 1, 10, 100, which exhibit a radial structure qualitatively similar 

to that of the corresponding fields in the plane C ■ 0. - The Figs. 25-27 

show E versus p for c ■ —0.9 end ut ■ 1, 10, 100; H ■ 1, 10, 100, 
z 

which is qualitatively similar co E^ in the plane c ■ 0. The Figs. 28- 
30 show versus p for C ■ -0.9 and ut ■ 1, 10, 100; H • 1, 10, 

100, which is qualiUtively similar to in the plane c - 0. - These 

comparisons indicate that the bulk discharge remains concentrated 
around the axis p ■ 0 from the cathode region C ■ -0.9 to the 
central region C • 0. It is remarkable that the discharge does not 
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spread significantly in radial direction with increasing axial distance 
C within the interval -0.9 < C S 0 although the ratio of anode and 

cathode radii is - 90. 

An cde Region, c ■ -K).9: In Figs. 31-39. ♦, E^. and V are 
represented versus p for C “ +0*9 ut “ 1. 10. 100; H ■ 1, 10. 

100. One recognizes that the discharge has spread radially, in 
particular at moderate values of wx and H. in the plane ? ■ 40.9 
due to the Influence of the anode at p • 0.9 and C » +1. - The 

Figs. 40-42 show versus p for C » +0.9 and oix « 1, 10. 100; 

H “ 1. 10. 100. For moderate values of ux and H. is strongest 
at p * 0,9. whereas E Is strongest at p « 0 for large values of 
UJX and/or H. - The Figs. 43-45 show versus p for i; - +0.9 and 
MX - 1. 10. 100; H - 1, 10. 100. Is most intr^nse i) at p « 0.9 

for moderate values of wx and H and 11) at p » 0 for large values 
of wx and/or H. - It is evident that the axial magnetic field 
inhibits the radial spreading of the discharge at sufficiently large 
values of ux and/or H. The discharge bends around toward the 
circular anode of radius ® 

anode plane C - +1» where Ac is the smaller the larger wx and/or 
H are. 

T!ie above results are readily applicable to ordinary dense 
discharges with known transport properties. An application of the 
theory presented to the thruster discharge is more difficult, since 
the "anomalous" transport properties (u»o.x) of low density discharge 
plasmas are not known. Other complications arise from the presence of 
the baffle electrode, which bisects the discharge region, and the 
inhomogeneity of the external magnetic field. With some confidence, 
however, the following qualitative conclusions should hold; 
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The torque produced by the Lorentz forces rotates the thruster 

_2 

discharge around the chamber axis. For |l| “1 amp» 10 Tesla • 

o-lO^ mho m"^, R - c - lO'^m, one has 0[v^] - lO'^m sec"^ and, hence. 

by Fig. 6, 0[V] - lO^m sec"^ for ux - 10, H - 10 . 

The usual assumption of an approximately homogeneous current 
distribution throughout the discharge space Is most probably Inapplicable. 
It must rather be assumed that the discharge current Is concentrated In a 
narrow region around the chamber axis and a thin layer on the cathode 
side of the baffle. In the annular gap surrounding the baffle, the 
current density should be extremely high and restricted to a thin layer 

adjacent to the baffle. 

Conductivity estimates based on the assumption of a homogeneous 
current density In the annular gap between the baffle and the cylinder 
piece will necessarily be too pessimistic. The radial spreading of the 
thruster discharge Is considerably reduced by the magnetic confinement 
field at Hall coefficients ojt»1. 
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FIG. 10: E versus p for C ■ 0, and wt ■ 1; 
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H » 1. 10, 100. 
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FIG. 11: E versus p for C ■ 0, and o>t 

z 

H - 1, 10, 100. 






PIG. 13: versus p for C ■ 0, and ut ■ 1 

H - 1, 10, 100. 
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FIG. 19: ♦, E^. and V versua p for C - "0.9 

and WT “ 10; H ■ 1. 
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PIG. 23: ♦, E^, and V versus p for C “ -0.9, 

and «T ■ 100; H ■ 10. 
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FIG. 26; E versu# p for ; - -0.9, and wx • 10 
z 

H - 1, 10, 100. 
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PIC. 44: versus p for C - +0.9, snd wt 

H • 1, 10. 100. 
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QUANTUM THEORY OP SPUTTERING: 
Application to Cross Section Detemlnstion 


Abstra ct 

Surface sputtering of polycrypt.illlne metals is explained 
theoretically by means of a 3-body sputtering mechanism involving the 
impinging ion and two metal atoms. By means of quantum-mechanical 
perturbation theory, a formula for the number S(E) of atoms sputtered 
on the average by an ion of energy E is derived from first principles. 
The theory agrees with experimental sputtering data in the low energy 
region above the threshold. As an application mercury-metal atom 
scattering cross sections are determined by quantitative comparison of 
the theoretical and experimental S(E)-values for sputtering mercury 
ions and various metals. 
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By means of quantxia^nechanlcal perturbation theory, the 
probability S(E) for an atom sputtered by an ion of energy E 
Incident on a polycrystalline metal is calculated for low ion energies, 
E j E^, where E^^ is the apparent sputtering threshold. Whereas the 
previous classical approaches to the analysis of sputtering 
contain phenomenological parameters (usually determined by fitting the 
experimented S(E)-curves) , the quantum-mechanical sputtering theory 
is based on first principles. The theory presented agrees with exper- 
imental sputtering data for low ion energies As an application, 

the formula derived for the sputtering ratio S(E) is used for the 
determination of the total scattering cross section for mercury atoms 

(recombined Hg-ions) interacting with atoms of various metallic solids 

, ^ , 11 , 12 ) 

from the corresponding experimental sputtering data. 

A binary collision between a surface atom of the solid and an ion 
Incident normal to the surface can evidently not lead to sputtering 
since the atom does not acquire a momentum component in the direction 
of the external normal of the surface. Similarly, sputtering is not 
likely to occur for smaller angles of ion Incidence if its energy is 
not large compared to the threshold energy for sputtering. It is evi- 
dent that sputtering, at ion energies of the order of the threshold 
energy, is a 3-body process involving one ion and two surface atoms of 
the solid. At higher ion energies, however, sputtering will result 
mainly from higher order many-body interactions. 

By restricting the theoretical considerations to ion energies E 
of the order of the threshold energy, E i E^. sputtering is regarded 

lon-atom-atom interaction. Furthermore, it is 


as the result of an 
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assumed that the solid is polycrystalline and has a sublimation energy 
which is on the average >, where the average is taken 

over the randomly distributed surfaces (ijk) of the crystallites. 

In this case, the sublimation energy E_ represents the average bind- 
ing energy of a surface atom of the polycrystalline solid. 

In the 3-body sputtering process, the incident ion transfers, on 
the average, the energy E (as well as kinetic energy) to the atom 
which is expelled and the energy o) 2E^ or B) 4E^ to the other atom 
depending on whether the latter is pushed to an a) unstable or 
6) stable interstitial lattice position. Accordingly, the threshold 
energies for the 3-body interactions a) and 6) are: 


^s 2E3 - 3 E^ , 


'6 


+ 4E^ - 5E^ 


Depending on whether the process a) or 6) occurs with dominant prob- 
ability, the apparent threshold (obtained by extrapolation of the 
experimental S(E) -curve, E -► E^) will be E^ ■ E^ or E^ Eg. If 
the cases a) and 3) have equal probability one might introduce an 
average threshold by 

% ^ 2 > - *^^8 • 


Indeed, some of the experimentally found thresholds E^(exper) 


11 , 12 ) 


can be explained by the theoretical formula E^ ■ ^^s* other cases, 
the formulae E - 3E and E. ■ 5E„ have to be used to explain the 

(X 8 g 8 

measured thresholds. This is demonstrated in Table I which com- 

pares the experimental“*~^ and theoretical threshold energies 

(E , Efl, E ) for different metals. Sputtering is in general not a 
Q 6 o 
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simple threshold process which can be defined by means of a single 
threshold value. This will be shown in detail through the following 

quantum-statistical considerations . 


TABLE I: Comparison of experimental and 

theoretical sputtering thresholds. 


Target 

Element 


Crystal 

Structure 


E (theor)[eV] 
o 


E^(exper) [eV] 


Ag 

Au 

Co 

Cu 

Fe 

Mo 

Nb 

Pt 

Ta 

Ti 

W 

Zr 


fee 

fee 

hep 

fee 

bcc 

bee 

bcc 

fee 

bcc 

hep 

bcc 

hep 


5^S- 

3E„ ■ 


16.75 

11.70 

22.00 

17.65 

20.60 

24.80 
30.84 

16.80 
> 24.00 
- 24.20 
• 35.20 
» 18.42 


17 

18 
22 
17 
20 

24 
32 
22 

25 
25 
35 
18 


In Tabl. I. the experimental threaholds have been taken from 
Stuart and Hehnerii’ . Theae author, concluded from their experi- 
mental data that the threahold la Independent from the mass 
ratio between the Incident Ion and target atom.ii^ and is. In 
first approximation, equal to - AEj. the average displacement 
threshold in radiation damage.ii^ It is seen that the agreeaent 
between the theoretical threahold. (E^.,.S„) *«<> the experiment, -1 
values E^(exper) la excellent, except in the ca.es Xu and W. 
Whether the n-procese or the 6-proce.s is dominant or both are 
(about) equally probable la apparently not dependent on the 
respective crystal system (fee, bcc, hep). 
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PERTURBATION THEORY 


In general, a sputtering Ion recombines with an electron into an 

9) 

atom as soon as it approaches the surface of a metal.~ This means 

that the incident ion Interacts actually like a neutral atom with the 

atoms of the solid. This neutralized ion is always referred to as 

**ion", in order to distinguish it from the "atoms" of the solid. 

Experiments Indicate that also the atom sputtered from the metal sur- 

9) 

face is electrically neutral.— 

When an ion of low energy as defined above hits the surface of 
a solid, one of the following processes may occur: 1) the ion is 

reflected without energy loss by the bound surface atom it encounters; 
2) the ion collides with a surface atom and quasi-simultaneously with 
a second atom so that 3-body sputtering results. The total probability 
for the ion to Interact in either of the two ways with the solid is 

Pjj - N^^^o(E) (1) 


where N is the number density of atoms in the solid and o(E) is 
the total (energy dependent) cross section for ion-atom scattering. 

Let Wj^(E) and W 2 (E) be the transition probability rates for the 
processes 1) and 2), respectively. The relative probability with which 
sputtering occurs is then 


W3<E) 


WzCE) W2(E) 

w^(E) + W 2 (E) * Wj^(E) ’ 


W2(E) << Wj^(E). 


( 2 ) 


Combining of Eqs. (1) and (2) yields for the sputtering rate, l.e., the 
number of atoms expelled on the average by one ion of energy E from 


the solid. 
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S(E) • o(E)N^^^ W„(E) . 

8 


( 3 ) 


On principle, o(E) can be calculated quantum mechanically, or 
determined experimentally. 

In the transition processes 1) or 2) , the ion Interacts with the 
surface of the solid within an area of the extension of the de Broglie 
wavelength, X « n//2mE. For this reason, the spatial part of the phase 
space is 

V - -y R^, R » R/(2mE)^^^ • (4) 

The transition probability w(E) from a state "i” to a state "f" is 
proportional to the matrix element |m^^ 1 in square and the density of 
final states dp/dE per unit energy,—^ 


w(E) 



2 ^ 
dE 


(5) 


where 


and 



ill 'I'f H d r 
0 


dE 


ft ,n d»(E) 


(2irR) 


dE 


( 6 ) 


(7) 


for a state containing n independent particles with moments 

p^, p^ p^ . 4»(E) is the volume of momentum space corresponding 

to the total energy E. H is the perturbation (operator) of the 
Hamiltonian of the ion-atom system which causes the transition i-^f, 
and i|f^ and are the wave functions of the total system before and 
after the transition which are normalized for the volume ft, ft > V. 
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Equations (5) through (7) represent the basis for the determination of 
the process probabilities Wj^(E) and w^CE)* 

1. Reflection State. Wj^(E) Is defined as the probability rate 

for the Ion to be reflected at the surface of the solid without energy 

loss. In the center of mass system, the Ion momentum Is p ■ *'^2mE In 

3 

the final state and the momentum space volume is tCE) * Airp /3. 
According to Eqs. (5) through (7), the transition probability for 
reflection is per unit time (n*l) 

Wi(E) [fi/(2TTR)^] (8) 

where is the matrix element of the transition 1). 

2. Sputtering State. w^CE) is defined as the probability rate 

for the 3-body sputtering state with threshold E^, o * a.S 

(E • 3E . E ■ 5E ) . In the center of mass system, the momenta of 
a s B 8 

the ion (i), the sputtered atom (s), and the second atom (a) can be 
chosen as 

Pi " ’^s “ "2 P ■ Pa 2 P ^ 

so that momentum is conserved ^jPj “ Since the potential energy 

E is expended in the sputtering interaction of type a - a, 6, the 
0 

total kinetic energy of the three particles is 

® * M ® 

Equation (10) represents an ellipsoid with the axes sections 
{4(mM/(m + 2M)]E*}^^^ and (ME*)^^^ in the slx-dlmenslonal space 
of the vectors p and q. Hence, the volume of the momentum space is 


73 


•<« - i 0 


( 11 ) 


From Eqs. (5) through (7) and (11) one obtains the transition probability 
for the sputtering state with threshold per unit time (n-2) 

«°(E) - f [!J/(2wn>^l^ »<= - 

where is the matrix element of the transition 2) and 

H(E - E^) - 1 or 0 for E > E^ + 0 or E < - 0 (Heaviside) , 

With the assumption w°(E) << Wj^(E), one obtains from Eqs. (2), 

(8), and (12) for the relative sputtering probability the expression. 


n l"if 


( 2 ) ,2 


.3/2 « - 


ft »< it / « 


H(E - E ) . (13) 

0 


iMif . 

The matrix elements in square are proportional to the 

probabilities for finding the interacting particles in the processes 1) 
and 2) in the interaction volume V (Eq. (1)], i.e., 


“ (V/£i)^. 


M 


( 2), 2 

if ' 


(V/fi)‘ 


(lA) 


since these are one and two independent particles in the interactions 
1) and 2), respectively. Substitution of Eqs. (14) and (4) into 
Eq. (13) leads to the following equation for the relative probability 


for sputtering with threshold E^ (o ■ o, 8): 

, JJl I (M/m) . 3/2 0 

24 4 + 2(M/m)^ -2 


H(E - E ) 
o 


(15) 


^2/1 ^ 1^/ 1^ Is s correction factor of the magnitude 

0 [h 2 /iJ ■ 1» which can be determined more accurately by evaluating the 
matrix elements if the force potentials of the interactions 1) and 2) 


are known. 
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SPUTTERING RATIO 


In the sputtering of a surface atom by an Ion, two fundamental 

cases (a) and (S) are distinguished which have the thresholds 

E - 3E and E ■ 5E , respectively. Let the probabilities for the 
as E s 

occurrence of the thresholds E^ and E^ be g^ and g^, which are 
normalized in the usual way. 




1 , 


ft ^ 
a , p 


. (16) 


It follows for the relative probability that sputtering occurs with 
either of the thresholds E^ and E^, 


W 


s 



(17) 


Substitution of Eqs. (17) and (15) into Eq. (3) yields for the number 
of atoms sputtered on the average by an ion of energy E: 


S(E) - o(E)N^^^ 


'1 + 


(M/m) ,3/2 

2(M/mV 


I \ 

o*a,B 


(E - E^)- 


H(E - 


V* 


(18) 


For applications, it is suitable to further simplify Eq. (18), 
which is strictly v.'*Ud only for ion energies E > E^ For 

example, if only one threshold E^ important 

(g « 1 or gg « 1) and the total scattering cross section o(E) 
varies slowly at E ■ (absence of resonances), Eq. (18) can be 


reduced to 


S(E, . o(Ey 


(E - EY 


H(E - E ) (19) 

o 


with 


0 (E) /E^ s o(E^)/E^^, t Z E^ . 


(20) 
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Equation (19) Is exactly of the form of the sputtering relation 

S(E) ■ const (E - E^) found phenomonologically by fitting experimental 

sputtering data.—^ It is commonly used in the extrapolation E -► E^ 

of experimental data to find the threshold Recently, also a 

relation S(E) - const (E - E^)^ has been employed in the extrapolatory 

12 ) 

determination of the threshold E , — which appears to be difficult 

o 

to justify theoretically. As expected, the thresholds reported in 
Refs. 11 and 12 are in general somewhat different. 
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APPLICATION 

In Che literature, measurements of the total scattering cross 

sections for Hg-atoms (recombined Ion) and target atoms such as Ag, Au, 

Co, Cu, Fe, Mo, Nb, Pt, Ta, Tl, W, Zr have apparently not been 

reported.—^ Theoretical cross section values are not available yet 

owing to the mathematical difficulties associated with the application 

, 1^) 

of quantum mechanical scattering theory to many-electron atoms. — 

For these reasons, the cross sections under consideration shall be 
determined here by comparing the theoretical (Eq. (19)] and experimental 
sputtering ratios S(E). 

In Figs. 1-12, the dashed curves repre.^enC the experimental 

12 ^ 

sputtering data of Askerov and Sena — for Hg ions and the (polycrlstalllne) 

tavget materials Ag, Au, Co, Cu, Fe, Mo, Nb, Pt, Ta, Ti, W and Zr 

(with the lowest S(E) value measured at E ■ E indicated by a dot). 

The corresponding theoretical sputtering curves S(E) , based on Eq. (19) 

and the theoretical thresholds given in Table I, are shown by solid 

lines. The cross section values chosen in such a way that 

the experimental and theoretical sputtering curves agree in the low 

energy region E <; fi, since theory and experiment should agree the better 

the lower the ion energy is (3-body sputtering model). The mass of Hg 

is m ■ 200.59 a.m.u., and the remaining constants M and N in 

Eq. (19) are given in Table I. The latter shows also the details of 

the calculation of the cross sections the experimental 

sputtering data by means of Eq. (19). It is seen that the cross 

sections o(E ) for atom-atom scatterings are between 10® to 10^ barns 
o 

at low energies, l.e. are of the order-of-magnltude expected (h 2 /j^ ■ !)• 



TABLE II 


Constants of Sputtering Formula and Cross 


Sections a(E ) for Various Target Atoms 
o 


Target 

AtortS 


Ag 

Au 

Co 

Cu 

Fe 

Mo 

Nb 

Pt 

Ta 

Tl 

W 

Zr 


Mlgr] 


. (M/m)^ .V^ 
4^-2M/m ^ 


1.7906 10 
3.2697 X 10 
9.7829 X 10 
1.0549 X 10 
9.2706 X 10 


-22 

I 

-22 

-23 


-22 


1.5926 X 10 
1.5422 X 10 
3.2385 X 10 


-22 

I 

-22 

I 

-22 


3.0037 X 10 
7.9514 X 10 
3.0519 X 10 


-22 

-22 


1.5143 X 10 


-22 


5.201 X 10 

1.856 X 10 

1.268 X 10 


-2 

I 

-1 

-2 


1.523 X 10 

1.111 X 10 

3.998 X 10 
3.716 X 10 


-2 

-2 

-2 

I 

-2 


1.820 X 10 


1.563 X 10 

7.581 X 10 


-1 

-1 

I 

-3 


1.615 X 10 
3.564 X 10 


-1 

I 

-2 




5.859 X 

5.903 X 

8.903 X 
8.468 X 
8,478 X 
5.657 X 
5.187 X 
6.599 X 

5.526 X 

5.659 X 
6.324 X 
4.253 X 


10 

10 

10 

10 

10 


22 

22 

I 

22 

22 

22 


10 

10 

10 

10 


22 


22 

.22 


10 

10 

10 


22 

I 

22 

22 




1.509 X 10 
1.516 X 10 
1.994 X 10 
1.928 X 10 
1.930 X 10 


15 

15 

15 

15 


1.474 X 10 
1.391 X 10 
1.633 X 10 
1.451 X 10 
1.474 X 10 


15 

15 

15 

15 

15 


1.587 X 10 

1.218 X 10 


15 

15 

15 


8.700 X 10 
2.672 X 10 


-2 


9.687 X 10 
1.993 X 10 
7.649 X 10 
2.248 X 10 


-2 

-2 

-1 

-2 

I 

-2 


3.719 X 10 
2.244 X 10 
9.593 X 10 


-2 

-2 

-3 


1.786 X 10 
1.344 X 10 
2.658 X 10 


-1 

-2 

-2 


1.384 X 10 
4.230 X 10 

1.166 X 10 

2.480 X 10 


-15 

-16 

I 

-15 

I 

-15 


9.512 X 10 

3.661 X 10 

6.417 X 10 
3.298 X 10 
1.587 X 10 
2.908 X 10 


,-16 

-16 

-16 

-16 


-16 


-15 


2.032 X 10 


5.236 X 10 


,-16 

-16 



I 
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The Figs. 1-12 indicate that the theoretical sputterlngs formula 

in Eq. (19) describes the experimental data rather well In the low 

energy region E ^ E^. The theoretical sputtering curves S(E) are 

plotted up to E - 120 eV, in order to show the deviations of Eq. (19) 

from the experimental data at larger ion energies, E. The 3-body 

sputtering model and the sputtering formula derived from it evidently 

represent adequate approximations only up to energies E ■ 2E^ to 3E^. 

It should be noted that the theoretical sputtering curves 

[Eq. (19)1 Are very sensitive towards changes in the thresholds E^. 

It can be shown that an adequate agreement between the experimental 

and theoretical sputtering curves S(E) cannot be obtained by choosing 

theoretical thresholds E noticeably different from those in Table I 

o 

and varying the values of the cross sections o(E^). Experimental or 
theoretical cross section values are obviously necessary to demonstrate 
the success of the quantum mechanical sputtering theory presented in a 


rigorous way. 
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Fig. 2: Theoretical ( — 

sputtering rati 
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Fig. 9; Theoretical ( ) and experimental ( 

sputtering ratios S(E) for Ta. 




0 20 4 0 60 80 100 * 20 

e[.v1 


Fig. 10; Theoretical ( ) and experimental ( 

sputtering ratios S(E) for Tl. 
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DEPOSITION OF SPUTTERING PRODUCTS ON SYSTEM SURFACES 


ABSTRACT 

An analytical theory is developed describing the deposition of 
sputtered atoms on system surfaces which cannot be seen along straight 
paths from the emitting surface. The boundary-value problem describing 
the diffusion of the sputtered atoms through the surrounding rarefied 
electron-ion plasma to the "hidden" system surfaces is formulated and 
treated analytically. It is shown that outer boundary-value problems 
of this type lead to a complex Integral equation, which requires 


numerical resolution. 
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In an ideal vacuum, sputtered atoms travel undeflected along 
straight paths determined by their initial velocities at the point of 
emission. Within this free particle flow, a system surface la reached 
by the sputtered atoms only if it can be seen along a straight line from 
the emitting surface. In reality, ion propulsion systems are surrounded 
by a very rarefied plasma consisting of escaped beam ions, recombined 
ions, and electrons. For this reason, always some of the sputtered 
atoms will be deflected out of their initial paths by Interacting through 
long-range forces (polarization forces) with the plasma particles so that 
they can reach system surfaces which are not seen along a straight line 
from the emitter. 

Figure 1 depicts the geometry of an idealized propulsion system 
which exhibits an emitting plane z“0, (accelerating grid), 

the rocket surfaces r "a, -c^z^O and z*-c, O^r^f^a, and the 
plane z • -d, a £ r i b of the solar energy collectors. All these 
system surfaces can be reached by the atoms sputtered from the emitter 
by diffusion through the rarefied plasma. The diffusion coefficient 
D is determined by the Vlasov equatlon^^ for the sputtered atoms inter- 
acting through weak long-range forces^^ with the plasma particles. In 
view of the mathematical difficulties aasoclated with the solution of 
outer boundary-value problems for the geometry in Fig. 1, a somewhat 
simpler system is studied here concxstlng of an emitting plane 
(z ■ 0, 0 < r a) , the upper rocket surface (r - a, -c £ z 0) and 
the pUne (z - -c, a < r < -) of the Ur energy collectors (Fig. 2). 
The latter is assumed to have infinite radial extension, r^^ ■ b • 



Plows of 
Osposition 
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.lac. 1« b». .ml b>>c.d (Mg. 1). Within th. «d.l of 
Fig. 2, p.ttlcl. d.p 08 ltlon on .p.tan .utfacs In the space a 1 -c 
cannot be analyzed. 


(MI 
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BOUNDARY-VALUE PROBLEM 

In the space z >-c, let the density of the sputtered atoms be 
designated by N(r,z) [cm“^J and the flux of emitted atoms at the 
emitter surface by I(r) [cm ^ • cm sec In steady state, the 

spatial distribution N - N(r,z) of sputtered atoms is determined by 
the boundary-value problem for the Laplace diffusion equation (Fig. 2): 

i!l! + 1 M . 0 U) 

.2 ^ r 3 t ,2 


where 


-1 


[3N(r,z)/3z]^_Q » -I(r)D H(a-r) 


N(r,z) - 0 , -c <z < 0 
r*a 


N(r,z) • 0 , a < r £ 
* z* -c 


. ( 2 ) 
, (3) 

. ( 4 ) 


and 

N(r,z) -► 0 , (r^ + z^) • » 

are the proper and improper boundary conditions, respectively. D 

designates the diffusion coefficient of the sputtered atoms in the 
rarefied plasma which represents a spatial average, D = < D(r,z) >. 
The Heaviside function is defined as 

H(a-r) - 1 , 0 < r < a * 

■ 0 , a < r i • 

The boundary conditions (3)- (4) imply that sputtered atoms 
arriving at the system surfaces are deposited there, l.e., do not 
return into the diffusion space. This assumption is at least 
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approxinately correct for nonheated surfacaa aa long ae the nu^er of 

atomic layers deposited Is not too large. The fluxes 

atoms arriving at the system surfaces r ■ a, -c <. z i 0 and z - -c. 


a < r < * are given by 


♦^(r - a,z) - -D 3N(r - a,t)/3r , -c 1 z < 0 

. (6) 

♦ (z - -c,r) - -D 3N(z - -c,r)/3z, a < r £ • 
z 

• (7) 

Accordingly, 


0 

ti ■ - 2iraD / [3N(r - a,z)/3r)dz 

-c 

, (8) 

00 

N - -2*D / [3N(r,z - -c)/3i]rdr 

Z--C , 

* (9) 


are the nuiAers of sputtered atoms deposited per unit time on the 
system surfaces r - a, -c <. z < 0 and z - -c, a ^ r < -, respectively. 
The above boundary-value problem can not be solved directly, l.e., 
requires a decomposition of the space z > -c Into appropriate sub- 
regions for which the associated boundary-value problems are readily 
solvable. In this approach the conmon boundary value at the decomposi- 
tion plane Is determined by an Integral equation. 


DECOKPOSITION BY ONE INTERFACE 
In Fig. 2, the space Is decomposed Into the regions 
1(0 <r<», 0_<z^») snd II (s ^ r ^ “C <. z 0). At the 

Interface, z ■ 0, a < r < •, the partUl 3M(r,z - 0)/az - ♦(r)D H(r-s) 
is Introduced as the conmon (unknown) boundary value #(r) of the 
adjacent regions 1 and II. 


99 


Region I . Nj(r,z) is described by Eq. (1)» and the boundary 
conditions (5), and 

3Nj(r,z - 0)/3z - -I(r)D"^ H(a - r) + i|<(r)D"^ H(r - a) 

Accordingly, a solution is sought In form of the Fourier integral 
(0 < z _< «) 


N (r.z) - 7 A(k)e"^ J^(kr)dk 
^ 0 


( 11 ) 


whence 

00 oo * 

- f k A(k) J (kr)dk - / J (kr) k dk |t-I(a)H(a - a) 

0 “ 0 ® 

(12 

+ ij;(a) H(a “ a)] J^(ka)ada 

by Eq. (10) in accordance with the Hankel transformation. Subatltutlon 
of the Fourier amplltudee A(k) from Eq. (8) Into Eq. (7) given the 
solution: 

N (r,z) - /[I(a)H(a - a) - iKa)H(a - a)] odo x 

^ 0 

7 J (kr) J (ko)dk * 

0 ® ° 

Nj(r,z) evidently satisfies Eq. (5) since J^Ckr) -►0 for r 
as well as Eq. (10) since 

CO 

( J (kr) J (ko) dk - 6(r - a)/o 
0 ® ° 

Region II . Nii(r.a) is described by Eqs. (1) and the boundary 
conditions (3)- (5), and 


I 


1 


I 
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aNjj(r.t - 0)/3« - 4»(r)/D ,« 1 r < - 

Accordingly ( a solution is sought in fom of the Fourier series 
(-C £ 2 i 0) : 


Nii(r,2) 


J B J (X r)e'^vC 
t- MOV 


V"1 


(14) 


(15) 


where 

J (X a) - 0 • 

0 V 

determines the eigenvalues v ■ 1,2,3»..> Equations (14)-(15) and the 
orthogonality relations for J^,(X^r) relate the Fourier coefficients. 

n - - . (17) 

to ^(r). Equation (15) evidently satisfies the boundary conditions 
(3)-(5). 

Integral Equation for <>(r) . In the solutions for N^(r.z) and 
Njj(r, 2 ), the function <»(r) is still unknown (Eqs, (13) and (15)]. 
Since the z-derivatives of Nj(r.z) and Njj(r.z) have already been 
»tched at the Interface, z - 0. a < r < •. ♦(r) is determined by the 
remaining continuity condition. 

Nj(r,z - 0) - Mjj(r,z - 0) . a < r < - • 

Subatituloa of Eqt.dS) »aa (15). (I7)yial4a for *(r) tha iotagral 
equation: 

m 

1 <»(o) H(o - a) Q(a,r) da • o(r) 

0 


( 19 ) 
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where 


Q(o,r) = a / J^Cha) 


«» tgh X c J (^^“) 

- I '' 2 ■ — ’ 


* v-l X^a Jl <>v®^ 


09 

o(r) = flM H(a - a)odo / J„(kr) J^(ka)dk , (2D 

0 0 

are the kernel and the souree of Eq. (19), respectively. Since 

j"j^(kr) J^(ko)dk - f [i K(4> "" - ^ K(%) H(o - r) . (22) 

where K(m - kh is the complete elliptic integral of the first kind. 
Q(a,r) and o(r) become: 


2 * tgh X c 

y ^ 2 a 2^ y — o — 

Q(o.r) - - - K( j) 2^ ^ ^^3 2 

r V A 1 V 


J (X a) , r ^ a , (23) 

o V 


,(r) - ^ fua) H(a - a) f K(^)da , r > a 
^0 ^ 


. (24) 


noting that o < a and r > a along the Interface of the regions I 
and II. and. hence, r > a in Eqs. (21)-(22). In case the sputtered 
ato« are emitted homogeneously, I(r) • 1„, E,. (24) reduce, to 

2 2 2 

2 “1 r»/® \ /I - ^ . r > a » (25) 

o(r) - - I r [^(-o) “ (1 - 2^ 2'^ * - 

' IT o y* r r 

where E(m - k^) is the complete elliptic integral of the second kind. 
Note that o(r ■ a) ■ 0 and o(r • “) • 0. 


I 
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From the mathematical point of view, Eqs. (13) and (15, 17), which 
give the atom density In the space a ^ -c, and Eq. (19), which deter- 
mines the remaining unknown function ^(r), represent the formal solution 
of the deposition problem. 

FURTHER REMARKS 

Attempts at solving the integral equation defined by Eqs. (19) 
and (23)- (24) in closed form were not successful. Consecutively, 
the boundary-value function iKz) was determined numerically from 
Eq. (19) by an iteration method. Substitution of a finite number of 
values ■ i|'(*j^)» 0 ^ < «, into Eq. (17) gave (approximate) 

numerical values for the Fourier expansion coefficients B^. The 
numerical solution M(r,z), obtained by substitution of these 8^ 
values Into the Fourier series In Eq. (15), however, showed Insufficient 
convergence, in particular at the comers of the system (Fig. 2). This 
difficulty is evidently due to (minor) numerical errors In the computation 

of B , which are added up In the Fourier series. 

V 

In order to find a simpler, analytical solution to the deposition 
problem, other analytical attempts were made, e.g. by decomposition of 
the diffusion space into three simple regions through two Interfaces 
(Fig. 3). However, this approach results In an Infinite system of 
algebraic equations and an even more complicated Integral equation. 

The analysis of the deposition of charged particles, produced by 
charge-exchange between sputtered atoms and beam Ions, leads to a 
similar outer boundary-value problem for coupled, nonlinear partial 
differential equations (collision-free electrohydrodynaalc equations 
and Poisson equation). For the latter reason, first the simpler problem 
of the deposition of neutral atoms by diffusion was ti^eated. 
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Figure 
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